function Fe = function_of_Fe(JB2)
%% 初始化 Fe
Fe = zeros(3, 1);

%% 当第 1 个节点处于第二类边界上时, 计算 Fe
if JB2(1,2) == 1
    kesi = -1; % 第 1 个节点的局部坐标
    fy = [1/2*kesi*(kesi-1);
        (1-kesi)*(1+kesi);
        1/2*kesi*(kesi+1)]; % 当 kesi = -1 时的插值函数
    Fe = fy*JB2(1,3); % Fe 的计算, 式(1.66)
end

%% 当第 3 个节点处于第二类边界上时, 计算 Fe
if JB2(1,2) == 3
    kesi = 1; % 第 3 个节点的局部坐标
    fy = [1/2*kesi*(kesi-1);
        (1-kesi)*(1+kesi);
        1/2*kesi*(kesi+1)];
    Fe = fy*JB2(1,3); % 式(1.66)
end

end